Polaronic effects in strongly coupled electron-phonon systems: 
Exact diagonalization results for the 2D Holstein t— J model 
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Ground-state and dynamical properties of the 2D Holstein t-J model are examined by means of 
direct Lanczos diagonalization, using a truncation method of the phononic Hilbert space. The single- 
hole spectral function shows the formation of a narrow hole-polaron band as the electron-phonon coupling 
increases, where the polaronic band collapse is favoured by strong Coulomb correlations. In the two-hole 
sector, the hole-hole correlations unambiguously indicate the existence of inter-site bipolaronic states. 
At quarter-filling, a polaronic superlattice is formed in the adiabatic strong-coupling regime. 



Polaronic features of dopant -induced charge car- 
riers have been detected in the copper-based 
high-T c compounds La2- x Sr x Cu04+ y , and even 
more in the isostructural nickel-based charge- 
transfer oxides La2- x Sr x Ni04+ y Q. To tackle the 
problem of (bi)polaron formation in such systems 
exhibiting besides a substantial electron-phonon 
(EP) coupling strong Coulomb interactions, it 
seems, at the moment, that approximation-free 
numerical quantum Monte-Carlo and exact di- 
agonalization (ED) analyses of generic model 
Hamiltonians yield the most reliable results. 
Along this line, by use of ED, the ground-state 
properties of Hubbard and t-J models with an 
on-site Holstein EP coupling have been studied 
on finite clusters in ID and 2D [2-4]. What is 
missing to date is an application of the power- 
ful ED technique to the calculation of dynamical 
properties of the Holstein t-J model (HtJM), in- 
cluding the full quantum nature of phonons. 

In this contribution, we employ the Lanczos 
algorithm in combination with a kernel polyno- 
mial moment expansion and the Maximum En- 
tropy method || to investigate the quasiparti- 
cle spectrum of a single hole-polaron in the 2D 
HtJM on a ten-site square lattice. Moreover, we 
compute different hole-hole/ phonon correlation 
functions at higher doping level in order to com- 
ment on hole-binding effects and charge-density- 
wave (CDW) formation. 

The HtJM is described by the Hamiltonian B 

H = H ph + H t _j-^/^J2( b i +b i)~ h i> W 



where TL p h and H t _ j represent the phonon part 
and standard t-J model, respectively, and the 
last term takes into account the interaction of 
doped holes (hi = 1 — X) ff cj^c^) with a single 
dispersionless phonon mode (which, e.g., may be 
thought of as representing a local apical-oxygen 



coupling; e p - EP coupling constant, 



bare 



phonon frequency). H acts in a projected Hilbert 
space without double occupancy. A general state 
of (Q) can be written as the direct product 1^) = 
J2i k c f IOeZ ® \k) P h, where I and k label the elec- 
tronic and bosonic basic states, respectively, and 

\k) P h = UlTiV^)- 1 ^ I0W Since the 
bosonic part of the Hilbert space is infinite di- 
mensional we use a truncation method 0] restrict- 
ing ourselves to phononic states with at most M 
phonons. To control our truncation procedure as 
a function of M, we calculate the weight of the 
m-phonon states in the ground state \\Pq) of H: 



„m\2 



Ez, fc Icf I 2 with m = E 



N 



In the 

numerical work convergence is achieved if the rel- 
ative error of Eq(M) is less than 10~ 7 . 

Figure g shows |c m | 2 for the 2D HtJM with a 
single hole at weak, intermediate and strong EP 
couplings [in what follows we have fixed J = 0.4 
(all energies are measured in units of t)]. The 
curves |c m | 2 (m) are bell-shaped and their max- 
ima correspond to the most probable number of 
phonon quanta in the ground state. These results, 
as well as the M-dependence of Eq at e p = 4 (see 
inset), confirm the importance of multi-phonon 
states in the (adiabatic) strong-coupling regime 
e p » 1, huj. 

In the analysis of the HtJM we start with 




Figure 1: Phonon- weight function \c m \ 2 and 
ground-state energy E (M) for the 2D HtJM. 

the study of just a single dynamic hole. Increas- 
ing the EP coupling in the adiabatic regime, we 
notice a continuous but rather sharp crossover 
from nearly-free polaron, described by an effec- 
tive transfer amplitude that is only weakly re- 
duced from its value in the pure t-J model, to a 
less mobile (small-size) adiabatic Holstein hole- 
polaron (AHP). Moreover, we found that the crit- 
ical EP coupling for the polaron transition is sub- 
stantially reduced due to prelocalization effects 
of the hole in the antiferromagnetic spin back- 
ground Q. To elucidate the difference between 
the FP and AHP limits and to demonstrate the 
formation of a hole-polaron band at large e p , in 
Fig. H we present the results for the i?-resolved 
spectral function 

A R {E) =Y.\^ N - 1) \^\< N) )\ 2 

n . a 

xS[E-(Eir i} -4 N) )]- (2) 

Of course in the very weak-coupling regime the 
spectral function is barely changed from that of 
the pure t-J model. Increasing e p , the lowest 
peaks in each Ag separate from the rest of the 
spectrum. These states become very close in en- 
ergy and a narrow well-separated energy band 
evolves in the strong-coupling case, where the 
gap to the next higher "band" is of the order of 
the phonon frequency fito. Note that the transi- 
tion to the AHP state is accompanied by a strong 
increase in the on-site hole-phonon correlations 
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Figure 2: Single hole spectral function Ag (E) for 
the 2D HtJM at Tiuj = 0.8 and e p = 0.5 (upper 
panel) 3.4 (lower panel), where the energy scale is 
shifted by Eq N ^ . The insets show the low-energy 
part of the spectra. The integrated density of 

states N(E) = J^dE 1 J2k A k( e ') ( bold solid 
line) is depicted in comparison with the result for 
the pure t-J model (dashed line). 

which are about one order in magnitude larger 
than the nearest-neighbour (NN) ones (cf. Fig. 
11 in Ref. (§). This indicates that the AHP 
quasiparticlc comprising a 'quasi-localized' hole 
and the phonon cloud is mainly confined to a 
single lattice site. As the phonon frequency is 
enlarged at fixed, the hole-phonon correlations 
are smeared out, and the crossover to the small 
hole-polaron is shifted to larger values of the EP 
coupling. 

Next we wish to discuss the two-hole prob- 
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Figure 3: Non-equivalent hole-hole pair corre- 
lation functions Cho-ho(\i — j\) in the two-hole 
ground state of the Ht JM as a function of e p ; here 
1-3 label NN, next NN, and third NN distances. 

lem. To get a feel for hole-binding effects, we 
have calculated the hole-hole correlation function 

C ho - ho (\i-j\) = (Ms P ,J)\kh\Me P ,J))- (3) 

Results for C ho _ ho (\i - j\) are given in Fig. |. 
In the weak-coupling region, Cho-ho{\i — j\) be- 
comes maximum at the largest distance of the 
ten-site lattice, while in the intermediate EP cou- 
pling regime the preference is on next NN pairs. 
As expected, increasing further e v , the maxi- 
mum in Cjuf—ho (\i - j\) is shifted to the short- 
est possible distance, indicating hole-hole attrac- 
tion. At e p ^> 1, the two holes become 'self- 
trapped' sharing a sizeable common lattice dis- 
tortion, i.e., a nearly immobile hole-bipolaron is 
formed. The behaviour of Cho-ho is found to 
be qualitatively similar for higher (lower) phonon 
frequencies (see inset), except that the crossings 
of different hole-hole correlation functions oc- 
cur at larger (smaller) values of e p , which shows 
the importance of both parameter ratios s p /t and 

Finally let us consider the quarter-filled band 
case. Here, we have investigated the more simple 
spinless fermion model (total S z — S'^ laa ,). In ac- 
cordance with previous approximative treatments 
based on the inhomogencous variational Lang- 
Firsov approach ||, we found, as the EP cou- 
pling increases, evidence for a transition from a 
FP state to a 2D polaronic superlattice, where the 
holes are self-trapped on every each other site. 




Figure 4: C ho _ ho (\6-j\) [left] and C ho _ ph (\6-j\) 
[right] are displayed at e p = 3 and Hu> = 0.8, 
where both diameter and gray level of the circles 
being proportional to the correlation strength. 

This crossover is signaled by a pronounced peak 
in the charge structure factor at (tt, 7f). To visu- 
alize the correlations in this state in more detail, 
in Fig. |] we have depicted Cho-ho{\i — j\) and 
the corresponding hole-phonon density correla- 
tion function C ho - ph (\i - j\) = (^olM^jl^o) as 
a function of \i— j\. Our exact results clearly show 
the phonon-dressing of the holes and the result- 
ing tendency towards CDW formation. A similar 
polaron ordering was observed in La1.5Sro.5MO4. 
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